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. ABSTRACT 

We present source counts at 6.7 /im and 15 /im from our maps of the 
. Hubble Deep Field region, reaching 38.6/iJy at 6.7 //m and 255/iJy at 15 

' These are the first ever extra-galactic number counts to be presented 

, at 6.7/Lim and are 3 decades fainter than IRAS at 12 /xm. Both source 

counts and a P{D) analysis suggest we have reached the ISO confusion 
limit at 15 fim: this will have important implications for future space 
missions. These data provide an excellent reference point for other ongoing 
ISO surveys. A no-evolution model at 15 /um is ruled out at > Sa, while 
two models which fit the steep IRAS 60 fj,m counts are acceptable. This 
provides important confirmation of the strong evolution seen in IRAS 
surveys. One of these models can then be ruled out from the 6.7 fim data. 

Key words: galaxies: formation - infrared: galaxies - surveys - galaxies: 
evolution - galaxies: star-burst - galaxies: Seyfert 



B. Mobasher^, C.P. Pearson^, M. Rowan-Robinson^, T.J. Sumner^, 
L. Danese^, D. Elbaz^, E. Egami^, M. Kontizas^, A. Lawrence^, 
R. McMahon^, H.U. Norgaard-Nielsen^, I. Perez- Fournon^^, I. Gonzalez-Serrano^^ 

^Astrophysics Group, Imperial College London, Blackett Laboratory, Prince Consort Road, London SW7 2BZ; 
^SISSA, Via Beirut 2-4, Trieste, Italy; 

^Service d'Astrophysique, Saclay, 91191, Gif-sur-Yvette, Cedex, France; 

Osservatorio Astronomico de Padova, Vicolo dell'Osservatorio 5, 1-35 122, Padova, Italy; 
^ Max-Planck-Institut fiir Extraterrestrische Physik, Giessenbachstrasse, D-8046, G arching bei Munchen, Germany; 
^Astronomical Institute, National Observatory of Athens, P.O.Box 200048, GR-118 10, Athens, Greece; 

Institute for Astronomy, University of Edinburgh, Blackford Hill, Edinburgh, EH9 3HJ; 
^Institute of Astronomy, The Observatories, Madingley Road, Cambridge, CBS OHA; 
^Danish Space Research Institute, Gl. Lundtoftevej 7, DK-2800 Lyngby, Copenhagen, Denmark; 

Instituto Astronomico de Canarias, Via Lactea, E-S8200 La Laguna, Tenerife, Canary Islands, Spain; 
Instituto de Fisica de Cantabria, Santander, Spain 



Accepted 1997 May 9; Received 1997 March 24; in original form 1996 Decmber 5 



1 INTRODUCTION 

IRAS galaxy surveys at 60 /im have consistently pro- 
vided good evidence for a population of star form- 
ing galaxies evolving with a strength comparable to 
AGN. This has been confirmed by numerous studies 
of count distributions and redshift surveys from 0.6Jy 
to 50mJy (Hacking and Houck 1987; Saunders et al. 



1990; Lonsdale et al. 1990; Oliver et al. 1995; Gre- 
gorich et al. 1995; Bertin, Dennefeld and Moshir et 
al. 1997). 

This evolving population discovered by IRAS 
could have very important implications for cosmo- 
logical studies. In particular these objects are likely 
to contribute strongly to the star formation history 
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of the Universe. [Other incidental issues include the 
possibly significant impact such objects could have on 
the cosmological far infrared background e.g. Oliver et 
al. (1992) and Franceschini et al. (1991).] The popula- 
tions seen by IRAS are mostly relatively low redshift 
{z < 0.2); deeper 750 surveys such as this provide a 
longer baseline in redshift giving a better handle on 
the nature of the evolution. 

This paper will discuss the source counts from 
our maps of the Hubble Deep Field (HDF, Williams 
et al. 1996). Our observations have been described by 
Serjeant et al. (1997; Paper I) and the source extrac- 
tion described by Goldschmidt et al. (1997; Paper II): 
a total of 27 sources were found at 6.7 /im, and 22 at 
15 fj,m. Further papers discuss the associations with 
optical galaxies (Mann et al. 1997; Paper IV) and the 
models for spectral energy distributions and the star 
formation history Rowan-Robinson et al. (1997; Pa- 
per V). 



2 OBSERVED SOURCE COUNTS AT 6.7 
AND 15 fiM 

In Paper II we detected 7 sources in our 6.7 fim maps 
and 19 in our 15 ^m maps using a well defined source 
detection algorithm. To convert these source lists into 
source counts requires an estimate of the area within 
which a source of observed fiux (S^) could have been 
detected. 

To estimate this we need to know the minimum 
flux (Slim) detectable at any position. The source de- 
tection algorithm requires m pixels to have intensity 
I{x, y) > T{x, y), where T{x, y) is the threshold inten- 
sity. The flux of resulting detections is estimated us- 
ing an estimated local background intensity B{x,y). 
Assuming a well determined point spread function 
(PSF), P{x — xo,y — yo), this algorithm gives us a 
detection limit 

Su^ixo, yo) = max. ( n-,y) - Bix,y) \ 
\P(x ~xo,y- yo) J 

where max. is a function giving the mth largest 
value. T is defined for various areas in Paper II and B 
is determined by running the sky annulus across the 
full image. 

The PSF is the only uncertainty in this estimate 
of the survey areas, we thus decided to investigate 
this in some detail. One estimate for P comes from 
the standard ESA products. The PSF is estimated at 
3x3 sub pixel positions for each filter and lens posi- 
tion. We drizzled these images together to produce a 
single PSF for each of our two observing modes. For 
both bands, the ESA PSF contains a large amount 
of flux in the wings. These wings may be caused by 
scattered light or other data reduction features in the 
ESA PSF. In any case we cannot use the PSF for 
analysing the effective area since only 20 per cent of 
the total flux is within the Airy disk, so we calculate a 



Table 1. Characteristics of various Point Spread Func- 
tions: column 1, peak intensity; column 2, full width half 
maximum; columnS flux within twice the nominal Airy 
disk diameter; column 4, intensity in mth brightest pixel 
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'revised' ESA PSF which is background subtracted in 
the same way as were the sources and normalised to 
1 within the source aperture. Since our observations 
involved long integrations, in which jitter might sig- 
nificantly blur the PSF, and also because our data re- 
duction did not take into account field distortions, we 
decided to estimate P empirically from the HDF ob- 
servations themselves. P was estimated by summing 
the intensities from a number of sources in an square 
aperture 7 pixels to a side (i.e. 7 arcsec at 6.7 /im and 
21 arcsec at 15 lira). The relatively small aperture at 
6.7 /xm was necessary to avoid including more than 
one source. At 15 /im we excluded objects near the 
boundaries and the fainter sources leaving a sample 
of 6 sources, while at 6.7 fira we used all sources in 
the complete sample. We then normalised such that 
^ P = 1 over the aperture. Some parameters of the 
PSFs discussed are summarised in Table ^ 

Figures ^ and ^ show the effective areas of the 
surveys to sources of a given fiux [^^(Siini < S)]. No- 
tice that that the curves do not pass through the 
origin in Figure |2|, this is because the sources were 
detected using a fixed global threshold but the back- 
ground is estimated locally, hence a source could in 
principal be detected in a high background region but 
then assigned zero or even negative flux. 

The faintest flux limit that can provide useful 
information is defined by the smallest useable area. 
We choose the smallest useful area to be 200 beams, 
using Figures |l] and ^ this translates to 38.6 and 255 
/ijy over 1.6 and 6.3 arcmin^ respectively. These flux 
limits include 6 of the 7 6.7 fira sources and 17 of 
the 19 15 /im sources. Allowing a smaller area of 40 
beams would have suggested a flux limit of 30.4 and 
161 /iJy in which case the faintest 15 fim source would 
be excluded and the faintest 6.7 fim source would be 
at the flux limit. 

At this point we can examine whether our HDF 
data approach the confusion limit of ISO. The max- 
imum area available at 15 fim is around 18 arcmin^ 
thus the complete sample from Paper II of 19 sources 
is below the classical confusion limit (more than 1 
source every 40 beams) see Figure |^. As most of the 
complete 6.7/im sample of 7 objects have fluxes are 
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0.00' 0.05 0.10 0.15 0.20 0.25 
S(6.7 ^/,m)/mJy 

Figure 1. Area in which the 6.7 /im map is sensitive to sources at a given flux. The dashed hne is calculated using the 
re-normalised ESA PSF. On the right is the number of sources that would cause one to exceed the classical confusion limit 
of 40 beams per source. 




0.0 0.2 0.4 0.6 0.8 1.0 
S(15 /_im.) / in Jy 

Figure 2. Area in which the 15 fim map is sensitive to sources at a given flux. The dashed line is calculated using the 
re-normalised ESA PSF. On the right is the number of sources that would cause one to exceed the classical confusion limit 
of 40 beams per source. 



above 40/.iJy, Figure |l| shows that they are confined to 
an area > 2 arcmin^ and thus are above the classical 
confusion limit. Including the supplementary list from 
Paper II provides a combined sample of 27 sources 
which could well represent a confused sample as the 
majcimum survey area available is ~ 6 arcmin^. 

We use the areas determined with the empiri- 
cal PSF to construct the observed differential source 
counts or integral counts 



j:dS 

j:>S 

The integral counts are shown in Figures |^ 
through 1^ together with models and (in Figures ^ and 
^ IRAS counts, which are discussed later. The inte- 
gral counts are shown rather than differential counts 
as these are more useful in practice, however the data 
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points are not independent and so Poisson la lim- 
its are shown as a hatched area. Rigorous statistical 
comparisons between data and models are best made 
by assessing the number of objects expected at given 
flux limits and for this we apply the area as a func- 
tion of flux to the model; these results are discussed 



in Section 4.4 



With the exception of the analysis of simulated 
data sets we exclude objects in our complete list that 
are not associated with optical counterparts in Paper 
IV. This criterion removes one source at 6.7 /im and 
5 at 15 /xm. 



3 STELLAR SOURCE COUNTS 

Before moving on to discuss the galaxy counts a brief 
word needs to be said about the stellar counts. At 
6.7 ^m there are no obvious stellar candidates in 
the complete sample. This is no real surprise since 
the IfDF was selected to exclude bright stars. Ex- 
trapolations from the models of Franceschini et al. 
(1991) predict 0.38 star arcmin"^ for 5*6.7 > 40/x Jy 
at these latitudes, i.e. 1.9 stars if the HDF area was 
not biased against stars, reasonably consistent with 
our finding none. At 15 fim we see one stellar image 
in the flanking fleld (KOHDFS J123709.8-f 621239). 
The Franceschini et al. (1991) stellar model would 
predict 0.1 star arcmin"^ i.e. around one to two stars 
in this field which is quite consistent with our single 
detection. This stellar object has been exclude from 
the galaxy counts. 



4 COMPARISON WITH FAR INFRARED 
SOURCE COUNT MODELS 

4.1 IRAS Galaxy Counts 

We can use the 12 /im IRAS galaxy counts to esti- 
mate the bright counts at 15 fim using the Rush et 
al. (1993) sample. In order to estimate the area of 
this sample we had to construct a mask to take ac- 
count of the IRAS missing strip and regions of high 
source density that may have been excluded or un- 
derrepresented. To this end we applied the QMW 
IRAS Galaxy Catalogue mask (Rowan-Robinson et 
al. 1991); this cut excluded only 43 galaxies, reducing 
the sample to 850 objects and provided an estimate of 
the area of 6.76 sr. These counts are shown in Figures 
I and |. 

4.2 Model Galaxy Populations 

We consider two models of galaxy number counts. 
Both models have well defined dust emission spec- 
tra specifically to accurately predict mid to far in- 
frared galaxy distributions. In addition both models 
include strongly evolving components sufficient to ex- 
plain the steep number counts at 60 fim. The specific 



populations and spectral energy distributions in the 
two models are however significantly different. Pear- 
son and Rowan- Robinson (1996; PRR) have described 
a galaxy population model involving five populations: 
normal galaxies; star-burst galaxies; hyper-luminous 
galaxies; Seyfert 1 galaxies and Seyfert 2 galaxies. We 
have predicted the counts at 6.7 ^m, 12 /^m and 15 nm 
using these models ignoring the hyper-luminous pop- 
ulation which will have negligible contribution. Here 
star-bursts and normal galaxies have 60 fim luminos- 
ity functions taken from Saunders et al. (1990), the 
Seyferts 12 fim luminosity functions come from Rush 
et al. (1993). Both Seyferts and star bursts evolve as 
L{z) = L(0)(l-|-z)^ The SEDs used for these galax- 
ies are based mainly on IRAS data and are described 
in Pearson and Rowan- Robinson (1996). The model 
assumes an f2 = 1 cosmology. These models were 
shown to provide a good fit to the /7?^5'60 /im counts. 
As they stand these models would be unable to ac- 
count for optical or K-band counts but would require 
a low fl or more strongly evolving star-burst popula- 
tion. The integral counts predicted by this model are 
shown in Figures ^ and ^ 

A second model comes from Franceschini et 
al. (1994; AF). The total counts are modelled as 
the sum of 5 populations: AGN; star-burst galaxies; 
spiral/irregular galaxies; SO galaxies; and elliptical 
galaxies. The late-type systems (Spirals, Irr and star- 
bursts) evolve as: L{z) = L{0)e^'^^^^ (where r is the 
look-back time), in an open universe (go = 0.15). The 
early-type systems (Ellipticals, SO) evolving accord- 
ing to Franceschini et al. (1994), i.e. assuming that a 
bright phase of active star-formation at 2 ~ 2 — 4 is 
obscured by dust quickly produced by the first stel- 
lar generations. This same model accounts in some 
detail for the sub-mm background as estimated by 
Puget et al. (1996). The AGN number density is set 
by the local IRAS samples at 12 fim (see Rush et al. 
1993). The evolution is calibrated so as to produce 
the hard X-ray background with a suitable distribu- 
tion of the dust /gas absorbing column densities. The 
integral countspredicted by this model are shown in 
Figures ^ and |^ 

The first impression from these Figures is of a 
surprising agreement between the models and the 
data, considering that in the case of the 15 /^m 
data the predictions were made from data 3 decades 
brighter, while at 6.7 /im there were no previo us d ata 



with which to normalise the models. In Section 4.4 we 
will compare the models and observed counts more 
rigorously, but first will discuss possible biases in this 
data. 



4.3 Simulations of the Observed Source 
Counts 

In any comparison of source counts with models, it is 
vital to discuss possible sources of bias that might be 
present. 
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Figure 3. 6.7 /im source counts from HDF (hatched region). Models at 6.7 fim based on Pearson and Rowan-Robinson 
(1996): all components solid; normal galaxies dotted; evolving star-bursts dash-dot-dot-dot; Seyferts 1 and 2 dash and 
dot-dash. The depths probed by other, forthcoming ISO surveys are indicated (ELAIS: Oliver et al. 1997; CAM-D &; 
CAM-S: e.g. Elbaz 1997; Taniguchi: Taniguchi et al. 1997) 
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Figure 4. 15 /im source counts from HDF (hatched region) with IRAS 12 fim counts (thick line) at bright end, from 
this paper and Rush et al. 1993 {IRAS data shifted to 15 fim using cirrus spectrum). Models based on Pearson and 
Rowan-Robinson (1996): all components solid; normal galaxies dotted; evolving star-bursts dash-dot-dot-dot; Seyferts 1 
and 2 dash and dot-dash. The depths probed by other, forthcoming ISO surveys are indicated (ELAIS: Oliver et al. 1997; 
CAM-D & CAM-S: e.g. Elbaz 1997; Taniguchi: Taniguchi et al. 1997) 



The first issue is the reliability of the detec- 
tions. Spurious sources arise particularly through 
non-Gaussian noise features remaining after the data 
reduction. Using optical associations Paper IV esti- 
mates lower limits to the reliability for the complete 
samples of 71 and 68 per cent at 6.7 /im and 15 fim. 
Thus the reliability could be comparable to the Pois- 
son errors and could be important to this analysis. 



For the direct comparison with the models we have 
excluded the non-associated detections, thus, by this 
estimate of reliability our list is 100 per cent reli- 
able. This method may introduce incompleteness, if 
we have thrown out genuine sources where the true as- 
sociation has an erroneously low likelihood (as might 
happen if our positional errors are not well under- 
stood). An alternative approach is to include all de- 
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tections and make a full assessment of the reliability. 
This requires a good estimate of the noise distribu- 
tion. 

Confusion between faint objects of high source 
density relative to the beam is an important factor 
in these surveys. Co-addition of overlapping source 
profiles means observed fluxes may be an overestimate 
of true single-source fluxes. On the other hand faint 
sources will tend to be merged together so decreasing 
their numbers. We can only estimate this bias with 
knowledge of the source counts at fainter fluxes. It is 
thus more appropriate to incorporate confusion to the 
models than to try and correct the data. 

Flux errors will cause objects to be scattered 
from fainter fluxes above the flux limit and vice versa. 
Since there are more fainter sources these are prefer- 
entially scattered over the flux limit into our survey. 
This 'Eddington' bias (Eddington 1913), also depends 
sensitively on counts at faint magnitudes and so will 
be included in the models rather than in the data. 

The total integral counts predicted at both wave- 
lengths from the first of these models have been used 
to create 100 synthetic gala:xy images using the em- 
pirical PSF. These source images were added to the 
background maps created with random pixel positions 
(see Paper I). These background maps include un- 
correlated noise and sky backgrounds but have real 
sources (or structure) suppressed to a negligible level. 
Some sources of correlated noise may also be sup- 
pressed in these background maps and this is a slight 
limitation to the simulations. 

These resulting images have been passed through 
the same source detection algorithm as was applied 
in Paper II thereby simulating all the possible biases 
discussed above. The only remaining possible source 
of incompleteness is from uncertainties in our PSF 
which may warrant further investigation. 

The average number of sources detected from 
these simulations are listed in Tables ^ and ^. They 
do not appear to be significantly different from the 
model counts on which they are based. Although the 
simulations are only made for the first model this has 
similar difi^erential counts to the second at faint fluxes 
so we would expect similar behaviour. In neither band 
do we see any appreciable discrepancies between the 
models and the simulations over the flux ranges in- 
cluding our data and conclude that any biases are 
small or cancel each other out. 

4.4 Comparison of Models 

Using the area of the survey to a given flux limit (Fig- 
ures l^and ^ and the number count models of Sec- 
tion we can estimate the number of sources of any 
given type that we would expect in this survey, these 
numbers are summarised in Tables Both mod- 

els are consistent with the total number of associated 
galaxies. The simulation at 15 /im predicts that there 
should be 8.37 objects but our non-associated list in- 



Table 2. Expected number of sources in ISO HDF Survey 
above 38.6;iJy at 6.7 /^ni, taking into account the areal 
sensitivity in Figure hi using the PRR model. The number 
of associated sources observed above this flux limit was 
5. For the simulation we did not exclude the one non- 
associated source above our flux limit since the simulation 
accounts for spurious sources. The second column indicates 
the number of sources expected at fluxes brighter than our 
brightest source. The final column is the probability of the 
observed fluxes being drawn from the model distribution, 
estimated using a K-S test 



Component 


N 


N 


K-S 




{S > SS.e^fJy) 


(S > 65.7/iJy) 


/% 


Normal: 


3.13 


2.27 


0.44 


Star-bursts: 


1.27 


0.92 


0.48 


Seyfert 1: 


0.66 


0.44 


1.30 


Seyfert 2: 


0.60 


0.42 


0.59 


Total: 


5.67 


3.95 


0.53 


Simulation: 


4.37 


2.86 


0.73 



Table 3. Expected number of sources in ISO HDF Survey 
above 38.6/jJy at 6.7 fim, taking into account the areal 
sensitivity in Figure ^ using the AF model. The number 
of sources observed above this flux limit was 5. The second 
column indicates the number of sources expected at fluxes 
brighter than our brightest source. The final column is the 
probability of the observed fluxes being drawn from the 
model distribution, estimated using a K-S test 



Component TV N K-S 

(S > 38.6/iJy) (5 > 65.7AtJy) /% 

Spiral: 1.33 0.78 3.56 

EUiptical/SO: 2.10 0.98 15.78 

AGN: 0.34 0.26 0.21 

Total: 3.77 2.03 6.87 



Table 4. Expected number of sources in ISO HDF Sur- 
vey above 255/jJy at 15 fim, taking into account the areal 
sensitivity in Figure |^ using the PRR. The number of as- 
sociated sources observed above this flux limit was 11. For 
the simulation we did not exclude the five non-associated 
source above our flux limit since the simulation accounts 
for spurious sources. The second column indicates the 
number of sources expected at fluxes brighter than our 
brightest source. The flnal column is the probability of the 
observed fluxes being drawn from the model distribution, 
estimated using a K-S test 



Component 


N 


N 


K-S 




{S > 255AjJy) 


{S > 726/iJy) 


/% 


Normal 


2.71 


0.94 


4.22 


Star-bursts 


2.81 


0.89 


5.83 


Seyfert 1 


0.73 


0.22 


9.00 


Seyfert 2 


0.78 


0.26 


4.57 


Total 


7.03 


2.32 


5.26 


Simulation 


8.37 


1.88 


2.31 
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Figure 5. 6.7 fim source counts from HDF (central hatched regions). Models based on Franceschini et al. (1994): all 
components solid; spiral galaxies dotted; E+sO dot-dash; AGN-dash. Faint end upper limits (hatched) come from the 
P(D) analysis. The depths probed by other, forthcoming ISO surveys are indicated (ELAIS: Oliver et al. 1997; CAM-D 
& CAM-S: e.g. Elbaz 1997; Taniguchi: Taniguchi et al. 1997) 
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Figure 6. 15 fira source counts from HDF (central hatched regions) with IRAS 12 /im counts (thick line) at bright end, from 
this paper and Rush et al. 1993 {IRAS data shifted to 15 /xm using cirrus spectrum). Models based on Franceschini et al. 
(1994): all components solid; spiral galaxies dotted; star-bursts dash-dot-dot-dot; sO dot-dash; AGN-dash Elliptical galaxies- 
short dashes. Faint end constraints (hatched) come from the P{D) analysis. The depths probed by other, forthcoming ISO 
surveys are indicated (ELAIS: Oliver et al. 1997; CAM-D & CAM-S: e.g. Elbaz 1997) 



eludes 16 galaxies. This is a marginally significant ex- 
cess and may indicate that either our association re- 
quirements are too harsh or that our simulated noise 
is not realistic; improved understanding of the ISO 
data in the near future will clarify this. Interestingly 
the PRR model predicts a significant fraction of ob- 
jects (~ 70 per cent) of the objects should be at fluxes 
brighter than our brightest source. This discrepancy 



is also seen in the integral count plots (Figure 




where we can also see a apparent difference in slope 
between the model and the data. 

To test whether discrepancies in the count slope 
are significant we perform a K-S test to assess the 
probability of the observed fluxes being drawn from a 
count distribution with the same shape as the models 
(or single components of the models), these probabil- 
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Table 5. Expected number of sources in ISO HDF Sur- 
vey above 255/^Jy at_15 fim, taking into account the areal 
sensitivity in Figure Busing the AF model. The number of 
sources observed above this flux limit was 11. The second 
column indicates the number of sources expected at fluxes 
brighter than our brightest source. The final column is the 
probability of the observed fluxes being drawn from the 
model distribution, estimated using a K-S test 



Component 


N 


N 


K-S 




(S > 255AtJy) 


{S > 726AtJy) 


/% 


Spiral: 


4.55 


1.59 


3.57 


Star-burst: 


1.26 


0.52 


1.25 


Elliptical: 


0.03 


0.01 


47.21 


SO: 


0.77 


0.13 


56.45 


AGN: 


0.25 


0.07 


14.93 


Total: 


6.87 


2.33 


4.76 



ities are given in Tables^,^|^|^ These results suggest 
that the PRR model at 6.7 /irn is ruled out with more 
than 99 per cent significance (using either the simu- 
lations or raw models), this is mainly due to the ex- 
pected fraction of bright objects. This model cannot 
be ruled out with more than 95 per cent confidence at 
the other wavelength although the simulation implies 
a more significant discrepancy, this may again suggest 
problems with our simulations or association criteria. 

In this test the AF model can also not be rejected 
with any high level of confidence (more than around 
95 percent) at either wavelength. 

It thus appears that these data are insufficient 
to rule out either model at 15 fim. The PRR model 
can be ruled out at the shorter wavelength due to the 
high fraction of bright galaxies predicted. It may be 
that revisions to the K-corrections of the "normal" 
component using improved SEDs may be sufficient 
to improve this model. Possibly problems with the 
"normal" component SED are also suggested by the 
noticeable over-prediction of the IRAS 15 ^m counts 
in the model. 

A "no evolution" version of the PRR star-burst 
component would predict only 0.28 galaxies 15 /xm i.e. 
4.5 in total. Since we comfortably detect 11 galaxies 
such a model can be ruled out at the 3 a level on 
integral counts alone. 



5 P{D) ANALYSIS 

Since our maps are close to the ISO confusion limit it 
is sensible to examine the low level fluctuations in the 
maps on a statistical level. This allows us to investi- 
gate source count models below the flux level at which 
individual sources can be resolved. To this end we ex- 
plore the distribution of deviations in flux in a given 
aperture from the mean level, the P{D) distribution. 
This analysis avoids the need for any source detection 
algorithm. The fiuctuations in a map about the mean 
intensity consist of two components. The first is the 



noise which will have both positive and negative val- 
ues, and will have some distribution function which 
we shall assume is Gaussian. The second is true fluctu- 
ations. These might arise from sources, cosmological 
background or Galactic foreground. In all cases this 
contribution will always be positive and thus skew a 
symmetric noise distribution. In this case we assume 
that any major asymmetric component arises from 
extra-galactic sources. 

Firstly we select sky regions and construct a his- 
togram of flux deviations in square apertures. These 
histograms are fitted with both Gaussian distribu- 
tions functions and the expected distribution func- 
tions from the source count models. These model dis- 
tribution functions are calculated following Frances- 
chini et al. (1989) and using the AF model for the 
source counts discussed above (the differential counts 
of this model are similar to those from the PRR model 
below the source detection limit so the AF model re- 
sults will be similar for both, and the model provides 
a better fit to the count distribution at 6.7 fim). The 
differential count distribution is first convolved with a 
model PSF (in this case a Gaussian PSF is assumed) 
to give the response function to single sources in the 
selected aperture. The P(D) can then be calculated 
assuming the sources are distributed on the sky as a 
Poisson process. 

Figure |^ shows the deflection distribution [D 
in fi3y) obtained from a 72x72 matrix of pixels de- 
rived from the inner portion of the drizzled mo- 
saic of Paper I. We estimate a sky standard devi- 
ation in the inner map outside obvious sources of 
(T = 16 — 17/iJy/aperture. (dotted line in the Figure 
The aperture is 6x6 arcsec^ area, enclosing a disc 
with diameter equal to the FWHM of the theoretical 
PSF. The mean intensity was 0.4 mjy arcsec"^. 

The continuous thick line is the convolution of 
the Gaussian noise with a model P{D) based on the 
counts appearing in Figure ^. The convolved curve 
provides a very good fit to the data (reduced ~ !)• 
A simple Gaussian cannot fit the data, even if the 
width is increased to a = 18 — WfiJy aperture"^ with 
reduced > 1-5. Thus there is a clear positive signal 
in the P{D). This means that the background in an 
extremely deep ISO exposure at 15/im is structured. 
Such structure is entirely consistent with being due 
to an smooth extrapolation of the ISO source counts 
observed above the flux threshold, and conflrms the 
counts. 

This P{D) analysis allows us to further constrain 
the models. The shape of the AF model counts was 
flxed and the normalisation allowed to vary. These 
constraints are illustrated in Figure |^ 

A similar analysis at 6.7 fim demonstrated that a 
simple Gaussian noise model was sufficient to provide 
a good fit to the observed P{D). This allowed only 
an upper limit on the normalisation of the counts to 
be determined. 
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Figure 7. P{D) analysis at 15 (see text) 



6 DISCUSSION OF POPULATIONS 

So far we have not made any use of the fact that 
our survey was conducted in the Hubble Deep Field 
and Flanking Fields where there is exceptional photo- 
metric and spectroscopy information available. These 
data allows us to determine the populations detected 
by ISO. It is still instructive to compare observed and 
predicted populations as this may provide clues for 
improvements to the models. 
One 

of the 6.7 urn sources (/5'0HDF2 J123646.4+621406) 
selected at 6.7 fim has broad emission lines (Paper 
IV). This is compatible with the predicted numbers 
of AGN in both of the models we have considered. Of 
the remaining 4 associated objects only one is com- 
patible with a normal cirrus spectrum galaxy and has 
an elliptical morphology, the others are better fitted 
with star-burst spectra and have spiral morphologies 
(Paper V). 

Excluding AGN, the first of the models (PRR) 
suggests that at 6.7 jj,m our sample should have a 
majority (71 per cent) of normal spiral galaxies with 
around 28 per cent star-burst galaxies. The expected 
probability of getting the observed one or fewer cirrus 
galaxies in a sample of 4 is ~ 7 per cent. Although 
the statistics are very poor this casts further doubt 
on the validity of PRR model at this wavelength. The 
classification of SEDs in Paper V is based mainly on 
the optical/IR luminosities. Extra information comes 
from the expected redshift distribution of these nor- 
mal galaxies and star-bursts which is shown in Figure 
^ A K-S test cannot rule out the possibility that the 
observed redshifts are drawn from the model distri- 
bution even if the redshift of the broad lined object 
is included: the probability of the data being drawn 
from the model distribution is 0.18. 

Paper V does not attempt to fit SEDs of the type 
used for the second model. However, this model pre- 
dicts that the non-AGN 6.7 fj,m sources would be 61 
per cent elliptical or SO and 39 per cent spirals. This 
is reasonably consistent with the morphologies and 
SEDs above. 

Only three of the 15 fim selected sources are as- 
sociated within the HDF itself. Two of these have spi- 
ral morphologies and are fitted by star-burst spectra. 



while the third has elliptical morphology and a cir- 
rus spectrum. With such limited statistics this is rea- 
sonably compatible with both models. Including the 
Flanking Field areas we find in Paper IV that six of 
the 15 fim selected sources have associations and red- 
shifts (or photometric redshifts). These are shown in 
Figure ^ together with the redshift distribution from 
the PRR model. Excluding the lowest redshift object, 
(which was excluded from our count analysis as be- 
ing below the 255mJy fiux limit that was also applied 
to our model redshift distributions) we find that we 
cannot rule out these redshifts being drawn from this 
expected distribution (probability 0.41). This test is 
hampered by small number statistics but also assumes 
that the objects in the fianking fields with redshifts 
are not biased in any way. Clearly obtaining spectra 
for the fianking field 15 fj,m sources is a high priority. 



7 IMPLICATIONS FOR FUTURE 
SURVEYS AND MISSIONS 

The fact that we have reached the ISO confusion 
limit at 15 fims with a sensitivity of around 0.2 
mjy has important implications for other surveys 
and missions, in particular NASA's Small Explorer 
Mission, the Wide-Field Infrared Explorer (WIRE, 
http:/ /www. ipac.caltech.edu/wire/). WIREis due for 



launch in September 1998 and plans to survey hun- 
dreds of square degrees at 12 and 25 /xm. The 
WIRE strategy is divided into three parts: Part I 
- a moderate-depth survey (60 percent of survey 
time); Part II - a deep, confusion-limited survey 
(30 per cent of the survey time) and Part III - an 
ultra-deep, confusion-distribution measurement. The 
areal coverage and integration times will be designed 
to achieve these aims and so the best estimate of 
the confusion limit is required prior to planning. 
Currently the WIRE team estimate confusion lim- 
i ts of between 0.067 mJy and 0.15 mJy at 12 n m 



( tittp: / / www.ipac.caltech.edu/WIRE/sensitiv.html ) 
depending on the evolutionary models. WIRE has a 
30cm mirror with poorer resolution than ISO so our 
results suggest that the confusion limit for WIRE will 
be brighter than these estimates. 
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7. 

Figure 8. Expected redshift distribution at 6.7 /im witliin tlie ISO HDF areas (accounting for the area dependent flux 
limit) for tlie star-burst and cirrus component of the PRR model. Over-plotted at arbitrary y position are the redshifts 
(photometric or spectroscopic) for the 5 objects with reliable associations in Paper IV. The object at a redshift of 0.96 has 
broad lines in the optical spectrum. The probability that these data are drawn from the expected distribution estimated 
using the K-S test is 0.54 or 018 if the broad-lined object is included. 
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Figure 9. Expected redshift distribution at 15 [ixa within the ISO HDF areas (accounting for the area dependent flux 
limit) for the star-burst and cirrus component of the PRR model. Over-plotted at arbitrary y position are the redshifts for 
the 6 objects with reliable associations and spectroscopic or photometric redshifts in Paper IV, the lowest redshift object 
has the lowest flux and was excluded from our count analysis. The probability that these data are drawn from the expected 
distribution estimated using the K-S test is 0.41, excluding the lowest z object. 



8 CONCLUSIONS 

We have investigated the source counts of our ISO 
observations at 6.7 and 15 /jm of the Hubble Deep 
Field. These are the deepest available in the mid to 
far infrared wave-bands, 3 decades fainter than IRAS 
at 15 and the first ever presentation of extra- 



galactic source counts at 6.7 /im. We reach 50/iJy over 
5 arcmin^ at 6.7 /xm, at 15 pm we cover 5 arcmin^ to 
a depth 200/.tJy and around 15 arcmin^ to 400 /iJy. 
In both bands we cover smaller areas to deeper flux 
limits. These data will be very useful reference points 
for ongoing ISO surveys e.g ELAIS (see e.g. Oliver et 
al. 1997), the CAM Deep Survey (see e.g. Elbaz 1997) 
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and the survey of Taniguchi et al. (1997). Both the 
source counts and a P{D) analysis indicate that we 
have reached the confusion limit of ISO at 15 /xm. 
This will have important implications for the con- 
struction of mid-infrared surveys with future space 
missions such as WIRE. At 6.7 ^m we are close to 
the confusion limit but have not yet reached it. At 
15 /im a no-evolution model is ruled out at > Scr pro- 
viding important confirmation of the strong evolution 
seen in IRAS surveys at a much greater depth. The 
counts appear to be steeper than expected from one a 
priori galaxy evolution model at 6.7 ^ra although the 
significance of this is dependent on our assumptions 
about the calibration and noise properties of the data 
and the FIR SED of the model which will improve 
with time. The P{D) analysis provides constraints on 
number count models below the level at which we can 
extract reliable source lists. 

Further information on the ISOHDF project 
can be found on the ISORDF WWW pages: see 
http: / /artcmis. ph. ic.ac.uk/hdf/ . 
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Figure 10. P{D) analysis at 15 /im (see text) 
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